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1. Introduction 

Although the issue of the deconfinement phase transition in strongly interacting matter can be 
addressed with first principle calculation in lattice field fheory, if is of inferesf fo see whaf can be 
obfained regarding the chiral symmetry restoration aspect of the transition within effective models 
sharing the same global symmetry as QCD. For vanishing chemical potential universality classes of 
the transitions of strongly interacting matter were predicted by effective models for three comers of 
the niu^d — m^-plane: the exact SU(3) chiral limit (mu^d,s = 0), the two flavor chiral case (triu^d = 0, 
ifis = °°) and fhe pure gauge fheory {mu^d,s = °°)- As a consequence of fhe compafibilify of fhe 
conclusions fhere has fo be also a fricrifical poinf (TCP) on fhe muj = 0 line [p. Apart from fhese 
corners and ^ = 0 axis no universal argumenfs exisf and fhe order of fhe fransifion for differenl 
values of fhe chemical potentials (baryonic, isospin) and quark masses, as well as fhe shape of fhe 
fransifion line befween regions of crossover and firsl order fransifions are a mailer of calculation. 

If is expecfed fhal al vanishing baryonic chemical polenfial for low values of the masses there 
is a line of second order phase transitions as a boundary of the first order phase transition region. 
Increasing the value of /ig the second order line sweeps along a critical surface. Theoretically it 
is important to map the phase-diagram as a function of the parameters of QCD, and also there can 
be experimentally accessible interesting points, like the critical end point (CEP) occurring at finite 
/ifi, for physical values of the quark masses. 

Apart from helping the mapping out of the QCD phase transition, the effective models can 
give insights in the interplay between chiral and 6 ( 4 ( 1 ) symmetry restoration, and how meson and 
baryon properties change during the chiral transition. They can tell us what is the soft mode at the 
CEP, and by predicting specific signatures they might help finding ifs acfual locafion [^. They can 
also give informalion on fhe non-equilibrium dynamics near TCP or CEP ||3|, 

Effective models play also an imporfanf role whenever particle physics is applied fo cosmo¬ 
logical or asfrophysical sfudies. In cosmology, if is in imporfanf fo know whaf fhe cold dark mailer 
is made of. The relic densily of one of fhe candidafes, fhe weakly inleracfing massive parlicles 
(WIMPs), depends on fhe variation of fhe number of relalivislic degrees of freedom and fhe energy 
densily of fhe universe. If was shown in [^] fhal if WIMPs freeze oul in fhe range of fhe QCD Iran- 
silion lemperature fhe prediction of Iheir relic concenlralion depends on fhe QCD equation of slate, 
which below Tc is faken from an effeclive model. In aslrophysics fhe slruclure of compacl slars is 
actively invesligaled due fo increasingly available new dala fhal indireclly provides information on 
fhe slate of mailer in Iheir core. Al presenl effeclive models are fhe only fools for determining fhe 
equation of slale of slrongly inleracfing mailer al very low temperature and high density. 

In the present article a brief overview of some of the effective models will be given, and then 
some recent results on the issues enumerated in the Abstract will be summarised. 


2. Effective chiral models 

A common feature of the effective models is that they are constructed using the global chiral 
symmetry of the massless QCD and its spontaneous symmetry breaking pattern SUi{3) x SUr{3) —> 
SUv{3). Eor equal quark-masses the four-divergence of the vector current vanishes, while the 
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divergence of the axial vector current vanishes only in the massless case: 


= iq[M,X^]q = qi{mi-mj)Xfjqj, (2.1) 

dp [q'fysX^'q] = iq{X‘‘,M}ysq = qi{mi + mj)X^jY5qj, 

where for three flavors M = and there is an order parameter (OP), the chiral con- 

1 2 

densate (0|[e^,<775A*^]|0) = --(OI^IA",A^}< 7 | 0 ) = —-5ab{0\qq\0), which when vanishes indi¬ 
cates the restoration of chiral symmetry. 

The Lagrangian contains explicit symmetry breaking terms of two kinds. One is the quark 
mass term, or in the linear sigma model the term containing the external sources, which gives mass 
to the pseudo-Goldstone bosons. The other type is the instanton motivated’t Hooft determinant 
which breaks the Ua{^) symmetry. 

The effective models are built upon mesons and/or baryons with the lowest mass and are 
expected to be most relevant in case of a crossover transition when the mesons are well defined 
degrees of freedom even above the pseudo-critical temperature or in case of a second order phase 
transition where the change of the order parameter is directly accessible. 

The effective models differ essentially from QCD in that they miss the mechanism for confine¬ 
ment and incorporate only the degrees of freedom thought to be relevant for the chiral transition. 
Because of this latter feature it is very important to have a faithful, fine-tuned parametrisation of 
any effective model. I mention, however, that recently there have been studies [^, [Q] that included 
confinement in chiral models through an effective potential constructed with its (approximate) or¬ 
der parameter, the Polyakov loop, in an attempt to give account of the lattice result which shows 
that both chiral and deconfinement transitions happen at about the same temperature [^. 

The chiral perturbation theory-ChPT 

Chiral perturbation theory [||] describes exactly the low energy dynamics of the Goldstone 
particles collected in a hermitian traceless matrix <I>(x): 


^{x) 


^n^jVx + rjslVs 71+ K+ \ 

71^ -7r°/\/2-hT]8/\/6 

V K -xmlVeJ 


( 2 . 2 ) 


The most general Lagrangian invariant under SUl{3) x SUr{3) chiral symmetry is written down in 
terms of the exponential of <I>, U (<I>) = appears as a systematic expansion in momenta 

and quark masses: 


+ UM^) + 




(2.3) 


In the general Ut{3) x Ur{3) case [IB] this is more complicated because each term of (^!^) can 
be multiplied by some function of rjo = —f/logdetC /s/6 (in this case U = What 

is important for us from the point of view of the parametrisation of the linear sigma model in the 
— m/f-plane is that ChPT determines the and mR dependence of the decay constants for 
pions (fjc) and the kaons (/r) as well as of 


fn=f 


1 - 


(2.4) 
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fK=f 


1 - 


P V4 

.( 2 ) 


■ {l^Ti + Mt] + ) — 4m^L4 — 4m^(L5 + 2 L 4 ) 


(2.5) 


= 2m|-3v|^^^+2(2m| + m^)(3vf^+ ^ %vf\2m\ + m]^){Li + 3Lf) 


16 


2 \2 


(M?7 “ — 4wi^/Xj^ + — { 6 Lg — 3L5 + 8L7) (m^j — niff) 

+ —L(i{m‘^ — 2m\ + rrf^rriff) + —L-i{m^ + 2rri^)^ 


(2.6) 


2 2 

where := ^2 2 so-called chiral logarithms, determined at an appropriate scale 

(e.g. Mo = 4?!/) by the masses of the members of the pseudoscalar octet. / is the pion decay 
constant in the chiral limit. The low energy constants L; are determined by the values of the masses 
and decay constants taken in the physical point. 

The Nambu-Jona-Lasinio (NJL) model 

The Lagrangian of the chiral SU (3) symmetric NJL model is written in terms of the current 
quarks with mass matrix m = diag(m„,mrf,mi) 

£ [{qX^qf + {qiysX^q)^] +gD{det[^(l + 75)^] + det[^(l - 75)^]}. 

(2.7) 

The determinant term breaks the 6(4 (1) symmetry while the second term is invariant under chiral 
U(3) group transformations. For two flavors the determinant term can be rewritten as 

2{det[^(l + 75)^] +det[^(l -75)^]} = {{qq)^ - {qZqf - {q^q? + {qiY5^qf], 


and by taking equal couplings go = gs ^ G one arrives at the original form of the NJL Lagrangian 

0 : 

= <?(«^ - m)q + G{ {qq)^ + {qiysfq)^}. (2.8) 

This is invariant under axial transformations ^ > (l — jia ■ f 75 )^ since qq ^ qq — a ■ qiTj^q and 
qiTYsq —> qiTysqff aqq. In the mean-field approximation the square of the quark bilinear is lin¬ 
earised with the replacement {qq)^ —> 2{qq){qq) — (qq)^ and the vacuum expectation value of the 
quark bilinear, the order parameter, determines the constituent quark mass M for which the self- 
consistent gap-equation reads: 

r M 

M = mud-2G{qq) =mud + 2NcNfG [1-nq{T,id)-ng(T,ld)]. (2.9) 

7|p|<A y/p2-hM2 


The natural, mesonic degrees of freedom K, o are obtained through the procedure of bosonisation. 
Monitoring the expectation value {qq), one can study the restoration of chiral symmetry. For a 
review of recent developments in the NJL model see |^] and references therein. 

The linear sigma model (LaM) 

The Lagrangian of the linear sigma model with symmetry breaking terms is constructed in 
terms of aN/ x My complex matrix M = TaiOa + illa) which contains scalar (a^) as well as pseudo¬ 
scalar (tIq) fields [|I3|]. For Ny < 4 the Lagrangian reads: 


= \.r{d^M^d^M + ixIm^M) - (M^M) - f 2 pM^M)^ -g(detM -f detM^) -f tr//(M-fM^). 

( 2 . 10 ) 


4 











Chiral phase transition in effective models of QCD 


Zs. Szep 


For A/y = 3 the generators Ta are proportional to the Gell-Mann matrices. There are two independent 
quartic couplings f\ and / 2 . The explicit symmetry breaking term contains the external fields ha 
{H = Taha) which gives mass to the pseudo-Goldstone bosons. The role of the determinant term is 
to break the chiral U {Nf) symmetry down to SU {Nf). 

In the broken symmetry phase the scalars associated with the diagonal generators display 
vacuum expectation values. Considering only the case of degenerate u, d quarks one has for two, 
three, and four flavors one, two, or three condensates, called non-strange, strange and charmed 
condensates, respectively. 

It is important to stress that the determinant term enters with a relative opposite sign into the 
expressions of the tree-level masses of scalars and pseudo-scalars: 


[m Si^)]ab — Mo 2 ^ab ^abc^c] abed SNf,4 ^abcd]^c^d'! 

[m Uo)]ab — '^ab ff SNf,3 '^abc^c] ff^Nf,4^abcd\^c^d- (2-11) 

The various tensors appearing in ( p.llD are explicitly given in []T4[]. The parameters of the model 
are determined using the zero temperature mass spectra and decay constants. For the case Nf = 3 
this will be explained in Section ^ 

0(N) model 

For two flavors the 8 d.o.f. of the complex matrix M introduced previously form a reducible 
representation of SUl{2) x SUr{2). The irreducible representations are given by two 0(4) multi- 
plets {o,n) and (ao,ri), whose masses are split by the determinant term: = 2g, — 

= 2g — 2f\ ao. Assuming that the (ag,?]) multiplet is heavier (formally g ^ °°) one can neglect 
it, (integrate it out) and one arrives at the 0(4) model. An equivalent formulation is to use instead 
of general complex matrices unitary matrices with positive determinant M = al + i^ ■ T written 
with the help of sigma and pion fields [151. If one increases the number of pions and rescales the 
couplings in such a way that the energy density is proportional to the number of d.o.f. per site 
(~ N^) and the masses stay finite (~ A°), one obtains a form of the Lagrangian suitable for a large 
N approximation (see the meson part in the Lagrangian ( 2.12 ), below). 

The chiral constituent quark model 

This model extends the LaM by including some effective fermionic degrees of freedom, the 
constituent quarks. For a large N treatment SUi{2) x SUr{2) LaM model is extended to Nf flavors, 
that is the number of pions is increased and the appearance of an other quartic coupling is disre¬ 
garded. The constituent quark mass is generated by the scalar condensate. In the broken symmetry 
phase characterised by the vacuum expectation value <I>, the Lagrangian is parameterised in a way 
to be used in a large N approximation. Fermions couple with a Yukawa coupling to the mesons: 


■^CQM — — 


A<I,4 + i„j2^2 

N- 

X ^ 9 

—<I>^ -|- m^<I> -|- h 

24 2 


6 


1 


VN + -[id^^a)^ + {d^7i^f 




X 


6Vn 


<I>[a^ -l-a7r''7r"] - 




24N 


+ ¥ix) 


id^Y^ -mq-^ (a(x) -h i^/WfYsT‘^K'‘ 


0) 




( 2 . 12 ) 


where niq = g<I> and <I>(0) = 12. This parametrisation was obtained by requiring that the fermion 

mass stays finite as A —> oo. We can see that due to the large N treatment the fermion coupling 
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to the pions is enhanced over the sigma. The fermion contribution is of order 1 / s/N and thus it 
precedes in an 1 /N expansion the NLO meson contribution, which is of order 1 /N. 


3. Solving the effective models, thermodynamics 

Since the perturbation theory fails to account for non-trivial features of phase transitions, re¬ 
summation of perturbative series is needed. Such method is the optimised perturbation theory 


(OPT) proposed in [16]. It consists of introducing into the Lagrangian a temperature dependent 
mass with the change —/Tq ^ M^{T) and treating the difference between this and the original one 
as a finite counter-term Anp- = M^{T) -I- /Tq : 


•i^’mass = -^M2(r)Tr<I>^<I>-h^AMTr<I>^<I>. 


(3.1) 


The introduced mass-term is fixed fhrough fhe Principle of Minimal Sensifivify by requiring fhaf 
fhe one loop pion-mass defined af zero exfemal momenfum, calculafed using free-level masses in 
fhe propagators appearing in fhe self-energy, equals fhe free level mass for fhe pions: 


Mi(T) = ,G->(p = 0,T)| =-m2(r)|„-A™2-£ 


kl 



= ml{T) 


tree 


(3.2) 

This resulfs in a self-consisfenf gap equafion for fhe pions, since all fhe ofher free level masses can 
be expressed wifh fhe pion mass and fhe vacuum expecfafion values: 


a=a,7i 


^\{T) —+2{2f\+f 2 )x^ + ^fiy^ + 2gy + ^ c„.Itp(mo;;(r)). 


(3.3) 


The summafion is undersfood over fhe members of bofh fhe scalar and pseudoscalar mulfiplefs. 


In Oder to defermine fhe order of fhe phase fransifion fhe gap-equafion (3.3) has to be solved 
simulfaneously wifh fhe equations of sfafe which read as: 


-e„- }xlx + 2gxy + Afixy^ + 2{2fi +/ 2 )-^ ^ ^ Jita\vp{ma,{T)) =0, (3.4) 

-Ey - f.ly + gx^ + Afix^y + A{fi +/2)j^ ■/4a,Itp('Wa,(r)) =0, (3.5) 


where fhe weighfs and of fhe fadpoles and fhe isospin mulfiplicify factors 7,- are lisfed in 
Appendix C of [|I^. The renormalization was discussed in [l^j (see also [jT/l]). Wifh fhe definifion 


of fhe mass given in ( |3.2D only fadpole fype infegrals appear, and fhe only complicafion is due fo 
fhe mixing secfor where fhe eigensfafes are fhe rj and rj', and a and /o respectively. Because OPT 
does nof modify fhe form of fhe relafions upon which fhe Goldsfone fheorem relies, fhis mefhod 


preserve fhe Goldsfone fheorem for fhe pions. However, as explained in [ |18| ] resumming only one 
paramefer in fhis model wifh fwo order paramefers has fhe consequence fhaf fhe Goldsfone fheorem 
for kaons is nof safisfied af finife T. 

Anofher (approximafe) solufion of fhe linear sigma model was given using fhe Harfree approx¬ 
imation wifhin fhe CJT formalism [19], see also [^]. Using Euclidean mefric, we illusfrafe fhis 


mefhod for fhe case of fhe 0(N) model wifh an explicif symmefry breaking exfemal field, in fhe 
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large N approximation. For the chiral SU (N) symmetry there is an additional complication due to 
the mixing sectors. The effective potential parameterised with the vacuum expectation value <I> and 
the dressed pion propagator Gjt, includes all the diagrams that do not become disconnected upon 
cutting two lines (in the Hartree approximation only the pion “double scoop” diagram arises) 

A^- 1 “ 


V[^,G^] =N 


—- 

2 24 


+ 


Jp JP 


+ 


24iv [Jp 


G::{<^-,k) 


(3.6) 


The vacuum expectation value and the dressed propagator are determined from the two stationarity 


conditions 


5V[^,Gn] 5F[<I>,G 


5G^ 


= 0 , 


5<I> 


= 0. The first one gives a Dyson-Schwinger equation 


G„\^-,k) = ^„-\^-k) + 


4 f G„{^-,k)=k^ + m^ + ^^^ + ^ f G^i^-k), 

o Jp b b Jp 


and since the self energy is momentum independent, with the parametrisation Gjt{^',k) = 


(3.7) 

1 


1 


T O ^ O f/ f 

one arrives at a gap equation for the resummed pion mass: Mt = m H—H— / -i: -y. 

6 bJpk^+M^ 


U^+Ml 


The second condition gives the equation of state: 



X 

f ^ 1 

m + 

6 

+ 6 

Jp k^+M^_ 


^ = h. 


(3.8) 


We see from the gap equation ( ^^ and the equation of state (3.8) that in the large N approximation 
that is Goldstone’s theorem is satisfied. This is not the case in general because it is 
known that without the large N approximation Hartree approximation has two problems: it gives 
first order (instead of second order) phase transition in the chiral limit and Goldstone’s theorem is 
not fulfilled. 

The lasf mefhod for solving an effective model which we discuss now is fhe large N approxi¬ 
mation applied here to the chiral constituent quark model of Sec. We have solved this model in 
the leading order of the large N approximation by taking the fermions into account perturbatively at 
one loop order although there is an infinite series of diagrams containing fermions and giving 
contribution at order 1 / s/N of the large N expansion. 

Imposing a gap equation for the pions which performs the resummation of daisy-type dia¬ 
grams one obtains the temperature dependence of the order parameter by simultaneously solving 
the gap equation and the equation of state: 


gap equation: 
EoS: 


m 2 = (p = 0) = m2 + |<f2 + ^Ts{M)-4^TF{mg) 

m2 + I (<1)2 + Tb{M)) - (m. 


^eff(‘I’) = ‘I’ 


-h = 0 


Goldstone 

theorem 

h 


<!)■ 


m 2 = 


Tbj are the temperature-dependent pion and fermion tadpole integrals, respectively {niq = g<I>). 
The form of the a propagator can be inferred from the consistency condition —iG^^{p = 0) = 


G,\p)=p^---2<^^ 


1 - |/b(p,M) 


(3.9) 
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where Ib^f are the temperature-dependent pion and fermion one-loop “fish” integrals. Analytically 
continuing the sigma propagator to the second Riemann sheet the couplings were determined by 
requiring that the ratio of the real and imaginary parts of the sigma pole on this sheet is close to 
phenomenologically expected ratio of the mass to width of the sigma particle nialTa ~ 1, while 
keeping the mass of the sigma in the range G (400,800) MeV. The Yukawa coupling g is 
fixed by fhe consfifuenf quark mass (one third of the nucleon mass) and g = m^/<I>(0) = 6.72. 
Consistency of the T = 0 a-pole with phenomenological values of and requires A ~ 400. 

Once the couplings are fixed one refurns fo fhe solution af finife femperafure/densify. 


4. Restoration of chiral and axial symmetry 


We sfarf fhe discussion wifh fhe effecls of fhe resforafion of chiral symmefry as felt by the 
sigma excitation. 

In an effective theory the sigma particle is the quantum fluctuation of the amplitude of the 
chiral order parameter, hence any change in the ground state of the system is reflected upon its 
properties. Sigma couples strongly to the pions, measurable through the a —^2n decay. Because 
of this, at r = 0 sigma shows up as a broad resonance in the scalar-isoscalar channel. Hatsuda and 


Kunihiro conjectured [ |22[ ] that since the sigma mass decreases during chiral symmetry restoration, 
the available phase space for a ^ 271 squeezes, therefore one has a chance to see the sigma as a 
sharp resonance in matter. Observing the changes in the sigma one can infer features of the chiral 
symmetry restoration. 


5.5 

— Re Po/f,j for the complex solution / . 

— solutions on the real axis of the 2 sheey 

5 

— solution on the real axis of the 1 sheet / 

4.5 

/ 

4 

:_ T* 

3.5 


3 

2m,(T)/f, t / 

2.5 

T** 

2 



0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 


T/m,,(0) 

Figure 1: The temperature dependence if the real part of the sigma pole. 


Following with the increasing temperature the trajectory of the complex sigma pole on the 
second Riemann sheet obtained in the 0(N) model, with the help of the LO of the large N approx¬ 
imation, one can see that the real part of the sigma goes below the value of the two pion threshold 
at T** ^ 0.69mj[{0) = 96.6 MeV, which means the suppression of the a —>^271 decay channel [^]. 
This happens at a lower temperature than T* ^ l.07mjc{0), where a is converted into a stable par¬ 
ticle. So, for T < T** sigma is indeed a sharper resonance compared with the vacuum case. This 


scenario differs from the one in [ 24 ] where the pole that goes through the threshold is not connected 
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with the complex sigma pole at T = 0. By assuming certain analytic properties of the propagator, 
we have shown that there are universal features of the spectral function when a pole corresponding 
to a particle at high temperature approaches in the complex energy plane with the variation of the 


temperature the threshold position of its two-body decay [|25[]. 

It is interesting to note, that the value of the temperature where nia = Inijc is quite close to 
what was found within the composite operator formalism applied to the QCD: T** = O.QSTc ~ 98 


MeV [26|. Another interesting point is that within the large N formalism the mass of the sigma 


appears to be rather low (at NLO it is about 350 MeV, []27p). 

The restoration of the chiral and axial symmetries was studied recently by many authors, using 


both the linear sigma model [28, 29],[14] and NJL [30]. The method of these studies was to identify 
the chiral partners and observe the convergence of their masses. The chiral partners are (7r,cj), 
(T],ao) in the case of chiral SU(2) symmetry and (7r,ao), {K, k), (t],C7) in the case of chiral SU(3) 
symmetry. The results of these investigations can be summarised as follows. 

The scalar masses decrease at first but then start to increase during the transition while the 
pseudoscalar masses stay constant for a long while then increase monotonically with the tempera¬ 
ture during the transition. The masses grow more rapidly when V/ increases, and there is a tendency 
of the chiral partners towards degeneracy at some T > Tc. In the absence of 6 ( 4 ( 1 ) anomaly Tc de¬ 
creases, and the degeneracy is more rapidly reached. The melting of the strange condensate is 
much slower compared with the non-strange one. 

For the effect of the (/(I) anomaly it was found that in the chiral limit with(out) (/a(1) anomaly 
the transition temperature increases (decreases) with V/ in contrast to (in agreement with) the 
QCD lattice results which according to [^] give for Nf = 2(3) Tc ~ 173(154) MeV. This signals 
that the anomaly is rapidly decreasing near the transition temperature. Indeed, this decrease can 


be observed, since the finite temperature lattice result on topological susceptibility [31] can be 
converted into the temperature-dependence of the strength of the determinant term, g£){T) by fitting 
the lattice result with the explicit formula of the susceptibility calculated in NJL ]^^. In this case 
the determinant term decreases not only due to the very slow melting of the strange condensate, 
but also due to the decrease of goiT)- Hence 6 ( 4 ( 1 ) symmetry is restored and the degeneracy 
of the SU(3) chiral partners is completed. Without the restoration of the 6 ( 4 ( 1 ) symmetry the 
degeneracy among parity partners arises only at very high temperature at which the determinant 
term would eventually vanish due to the vanishing of the strange condensate. (We have to recall 
that the determinant enters with opposite sign into the expressions of the scalar and pseudoscalar 
masses.) 


5. fiB — T phase diagram for Nf = 2 

The left hand side figure of Fig. ^ shows a generic phase diagram typical for those obtained 
with effective models. In particular this figure was obtained for the chiral limit, solved in the 
LO of the large N approximation, with the fermions taken into account perturbatively at one 


loop order ]21]. The emerging diagram is qualitatively correct, but when compared with lat¬ 


tice results, one observes that the transition temperature is too low (the lattice result of [33] is 
Tc{flB = 0) = 173(8)MeV) and the present estimate for TCP is presumably located at a too high 
chemical potential. The transition line, which is a parabola, and the location of TCP can be obtained 
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Figure 2: L.h.s.: The jJ-s — T phase diagram of the constituent quark model in the chiral limit. 
Comparison of the CQM result from with lattice results (from top to bottom): Nf = 2 Nf = 
= 4 1361 (the original figure was taken from [ 351). 


R.h.s: 

3@, 


analytically in terms of polylogarithms: 


order+spinodal line 
order line ends when 


>n^+ ( \ 
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A 

6 471 ^ 




2n^Nf 

4“ Li« 


^. = 0 , 


, cTc 
— In — 
«=o Mq 


= 0 


Ttcp 

l^TCP- 


These formulae are obtained by expanding the equation of state in a power series of the expectation 
value and by taking the coefficient of the linear and cubic terms, which corre^ond to the quadratic 
and quartic terms of the effective potential. In the formulae above —z = e^ (the fugacity), c = 
Nc = 3, and Mq is the renormalisation scale. 

The curvature of the line of second order phase transitions at tto = 0 is = 0.101, 

which is somewhat higher than what is obtained on the lattice. (The masses are different, strictly 
speaking there is no much sense in comparing). This can be seen also in the r.h.s. of Fig. § taken 
from [|^] in which we inserted our phase transition line (the green line which ends in the lower 
right comer of the figure). The reason for this is that the large value of the Yukawa coupling lowers 
the critical temperature at vanishing chemical potential and bends too strongly the transition line. 
We have checked that dynamical generation of the fermion mass, e.g. solving a gap equation for it, 
lowers the effective Yukawa coupling. This result and the fact that according to Fig. ^ the line of 
the reduced temperatures lies within 10% of the lattice result, mean that it might be worth thinking 
about a more complete inclusion of fermions, e.g. beyond one loop order. 

By looking at different effective approaches published in the literature (see a list in [^7|]) we 
can see that all give low values for the temperature and high values for the chemical potential at 
TCP or CEP when compared with the result of Fodor and Katz ||^]. This may be partly due to the 
fact that these results are for two flavors and not 2+1 (it was emphasised in Ref. how sensitive 
is the dependence of the critical value of /Tg on the value of the quark masses). But this also raises 
a conceptual problem of whether the QCD phase transition is driven by the chiral symmetry. 
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As was shown recently in [^] using renormalisation group techniques the inclusion of gluon 
degrees of freedom proves to be important at T 150 MeV if one wants to compare the equation 
of state of the effective model with the one obtained on the lattice. 

The conclusion of the resonance gas model (which includes 1026 resonances) is that the phase 
transition is driven by the higher excited hadrons rather than by the light degrees of freedom [pT|]. 
This model managed to reproduce the temperature dependence of the energy density measured on 
the lattice. It is interesting that a pion gas would give only 15% of the critical energy density at the 
transition temperature. Imposing constant energy density the resonance gas model also reproduced 
the phase boundary measured on the lattice using Taylor expansion in /Tg. 

The discrepancy between the freezout curve and phase boundary obtained with the resonance 
gas model shows that there are strong interactions between baryons and mesons when the density 
of the baryons exceeds the density of mesons and without interactions this model is not able to 
account for how criticality occurs. 

I conclude that the investigation of the phase-diagram in the /Tg — T plane with an effective 
3-flavor model is definitely of interest. 


6. The nature of the soft mode 


In the physical point of the quark mass-plane a critical end point was found on lattice [38] at a 
finite value of the chemical potential and the question is what causes criticality. At the critical end 
point two minima and a maximum meet which results in a very flat effective potential of the chiral 
order parameter <I>. This is shown also by the diverging peak of the scalar susceptibility d^/dh. 
This flattening of the effective potential does not necessarily imply the appearance of a gapless 
sigma mode (sigma meson with vanishing mass). This is because the non-interchangeability of 
two limits ' : q-limit {co 0, then q ^ 0), which applies to static quantities as the curvature of the 
effective potential, and m-limit (q ^ 0, then co —>■ 0), which applies to the study of a sigma pole at 
rest on the second Riemann sheet (resonance). 


It was argued in Ref. [42] that what is important from the point of view of the signature 
of the genuinely singular end point is the dynamical feature rather than the static one. The time 
during which the correlation length ^ reaches its equilibrium diverges as T ~ (critical slowing 
down), where z is the dynamic scaling exponent. The finite evolution time limits the correlation 
length to (§ < (time)'/^. Because z > 1 and in heavy-ion collision the evolution times are of the 
same order as the spatial size, it is the finite time of the evolution rather than the spatial size 
limitation <size), which imposes a limit on the value of ^. This is why, it is important to find 
the dynamical universality class CEP belongs to, according to the classification of [|4^, and also 
the correct identification of the soft mode, the mode whose characteristic frequency vanishes at 
CEP, and which determines the characteristic time of the evolution. 

I refer the reader to the literature in questions which concern the hydrodynamical nature of the 
soft mode, in presence of conserved baryon number and energy densities, and on the demonstration 
that the dynamic universality class of CEP is that of model H, i.e. the liquid-gas phase transition 


'Using the relevant finite temperature expressions, one can easily prove that although, by definition, the bosonic 
one-loop bubble integral with equal masses is related at zero external four-momentum <? = (®, q) to the bosonic tad-pole 
integral through leiq = 0,M) = dTslM)/dM^, this relation holds only in the q-limit, not in the co-limit. 
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, 142], [04], Here I will discuss what can be found about the nature of the soft mode with model 


calculations focusing on the scalar order parameter. 

In the 0(4) model without explicit symmetry breaking there is a genuine second order tran¬ 
sition accompanied by a change in the symmetry of the ground state. Here, what is responsible 


for the critical behaviour and dynamical scaling is the pole in the sigma channel [|45[]. Following 
the trajectory of the sigma pole (corresponding to the particle a rest) of the analytically continued 
Green’s function on the second Riemann sheet, it turns out that the pole arrives on the imaginary 
axis, splits up into two poles, and the one that goes to the origin will produce critical scaling. In the 
vicinity of the critical point the dynamical scaling and the equation for the soft mode were given in 
p5| ] (see also the references therein). 

In the 2 flavor constituent quark model in the chiral limit, along the line of the chiral critical 
points, the soft mode is also the sigma meson, which becomes the degenerate chiral partner of the 
massless pions. The sigma pole on the 2nd Riemann sheet goes to the origin as the order parameter 
<I>(r) goes to zero. The real and imaginary parts of the pole scale as Repo <I>(r),Impo~<I>^(r), 
where <I> ~ (T — Tc)^. For (5 the mean-field exponent was obtained [1\] : jS = 1/2 for the second 
order line and jS = 1/4 at the tricritical point. The case of TCP is interesting because it prompts 
the question of what kind of singular behaviour might appear additionally there, where the sigma 
meson mode is already gapless (see below). 

For the physical pion mass and two flavors no symmetry change of the ground state accompa¬ 
nies the second order transition at CEP and in consequence from the point of view of the symmetry 
the existence of CEP is accidental. Both in CQM [ p^ ] and NJE [47] the sigma pole stays massive, 
and a peak develops in the scalar isoscalar spectral function in the space-like region (O < |q|. The 
location of its maximum goes to zero as the momentum goes to zero and it diverges in the q-limit. 

In the context of the NJE model what is responsible for the behaviour seen in the spectral 
functions is an imaginary pole of the scalar susceptibility iXmm = = d{qq) jdmud, 

where F{T,}x) is the free energy density) that goes to the origin as |q| — > 0 (diffusive mode) [El- 
In consequence the soft mode is the fluctuation of the scalar density. To study the relative weight 
of the diffusive mode in the susceptibility, the ratio R = (Zmm(0,0+) — Zmm(0^,0))/Zmm(0,0+) 
was introduced, where the numerator is the difference between the q-limit and m-limit. When the 
diffusive mode is absent R = 0. Because there is no contribution from the diffusive mode in the 
second term of the numerator, when the diffusive mode dominates R ^ 1. By studying this ratio R, 


it was shown in [44] that the switch of the soft mode from the sigma meson mode to the diffusive 
mode (hydrodynamical mode) occurs at the tricritical point, where R = 1. 


7. The boundary of the first order phase transition in the m^z — m/(:-plane 


m. 


There are several lattice studies of thermodynamics with varying degenerate quarks mu = md = 
which show that the value of ttip^ on the boundary drops substantially when improved lattice 


actions are used, namely from 290 MeV [^] or « 270 MeV p9|] to = 67(18) MeV 


[ 50 ] or ever further down to < 65 MeV [51] (see also the talk by O. Philipsen at this meeting). In 
view of this result we can think that the linear sigma model might complement lattice investigations 
if its coupling parameters can be determined accurately because the lighter are the pseudoscalars 
the better should become the approximate description based on EaM. This model was claimed to 
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give a low value for the critical pion mass: = 51 MeV [ ]52| ] or = 47 MeV [53] (see also 

[p^). In these investigations only the Goldstone mass was tuned by varying the explicit symmetry 
breaking terms and the other parameters were kept fixed at their values determined in the physical 
point. 


We have used ChPT when moving away from the physical point | |18| ] for the determination of 
the pion and kaon decay constants and +m^,, to obtain the value of the couplings at each 

point of the — m/^-plane (see Eqs. (p!^), (p!^, and (^)). Below I present schematically how 
the process of the parametrisation is performed. 


input: 


output: prediction: 



Ey=0 


The strange and non-strange condensates are determined through the PCAC relations x = f^, 
y = {IfK — f 11 ) 1 ^/^, while the masses determine the couplings (for the explicit formulae see [^] 
or l]!^). The predicted values of and are in remarkable agreement with the values of these 
masses given by ChPT up to = 800 MeV, even for vanishing pion mass. Unfortunately, at 
tree-level, in the pseudoscalar sector one of the quartic couplings, f\ and the mass parameter /Tq 
appear only in the combination = —/Tq +4/i (x^ +y^). To disentangle them we have to use also 
the scalar sector. Since the dependence of scalar masses on mjc and is not known, we need to 
make assumptions, which diminish the predictive power of our solution. This can be seen in Fig. || 
which shows the boundary of the first order transition region on the (m^ — mjf)-plane for various 
assumptions made for the continuation of the scalar mass spectrum. We expect that the exact phase 
boundary curve lies in the light-grey shaded region. 


Our rough estimate for the critical pion mass on the diagonal is merit (diag) = 40 ± 20 MeV, in 
nice agreement with the results of the latest effective model and lattice studies. In order to be able 
to use the tree-level parametrisation we used a quasi-particle approximation, that is we neglected 
the finite parts of the T = 0 quantum fluctuations. We note, that we were not able to locate the TCP 
on the = 0 axis. Recently a work performed in the NJL model found the TCP at the critical 
value of mj = 10 MeV [Q], which is much lower than expected when extrapolating lattice results. 
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Figure 3: Phase boundary for different alternatives. Al: it is assumed that there is no mixing in the scalar 
x — y sector (m^^ = 0) A2: it is assumed that the SUl{3) x SUr{3) Gell-Mann-Okubo mass formula for the 
scalars is satisfied in the full (m® — m/f)-plane with the same accuracy as at the physical point. We consider 
complete U(3) ChPT and the leading order large-Ac ChPT at Crossover (first order transition) takes 

places at the right (left) of the shaded region. 


8. Conclusions 


In this contribution we reviewed some of the results obtained with effective models of QCD 
on the chiral symmetry restoration. We have shown that one can parameterise the linear sigma 
model in the entire — mi^-plane employing the results of the chiral perturbation theory. We 
have determined the boundary of the first order phase transition line which is below = 65 MeV. 
For a more reliable determination of the boundary, information on the quark-mass dependence of 
the scalar sector is needed. 

We have shown that effective models correctly describe the changes of the meson properties 
during the chiral transition and indicate that 6 ( 4 ( 1 ) has to be restored as a prerequisite for restoration 
of SU (3) chiral symmetry. The sigma meson is the soft mode at the critical point of the 0(4) model 
with massless pions, and in both the 2 flavor constituent quark model and NIL model in the chiral 
limit along the line of the chiral critical points, but stays massive at CEP. At CEP the role of the 
chiral symmetry is rather secondary, here the soft mode has a hydrodynamical nature. 

An analytical determination of the fl — T phase boundary was given for two flavors in the chiral 
limit. The obtained phase diagram is qualitatively correct, but the transition temperature is low and 
the value of the chemical potential at TCP/CEP is apparently too high. Eight mesonic degrees of 
freedom, usually included into effective theories of QCD do not seem to be the single driving force 
behind the transition. In view of this, the extension of the 3-flavor chiral quark model to finite 
density is of imminent importance. It would be also interesting to develop large N techniques for 
SUl{N) X SUr{N) chiral models, since this approach has some advantage when used together with 
other resummation techniques (2PI or Dyson-Schwinger). 
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